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CONDITION ANALYSIS 



Background of the Invention 

The present invention relates to a method and apparatus for determining a program frr a subject, 
5 and to a method and apparatus for determining subject parameter values representing the effect of a 
condition on the subject. 

Description of the Prior Art 

The reference to any prior art in this specification is not, and should not be taken as, an 
10 acknowledgement or any form of suggestion that the prior art forms part of the common general 
knowledge. 

Currently, it is known to determine medication programs to allow drugs to be administered for 
different medical conditions. However the determination of the drug programs typically requires 
15 years of experimentation. Even then the regimes are typically fairly simplistic and rely on the 
patient taking specified quantities of medication at various time intervals- 

This approach to medication administration suffers from a number of drawbacks. Firsdy, the 
medication regimes are usually general to a condition and not specific to a respective patient. This is 
20 important as different patients often respond differently to the same medication regime. Tailored 
medication regimes are sometimes provided. However, this is currently only possible after 
monitoring the effect of the empirical medication doses on the patient, which can have a detrimental 
effect on the patient's health. Furthermore, not all effects of the medication are always readily 
apparent, so it is not always possible to determine exactly what the effect of the medication has 



This is exacerbated by the fact that assessment of the effect of medication is usually subjective, and 
is not quantified to allow the effect of medication to be determined absolutely. 

30 Secondly, it is typical for only a limited number of factors to be taken into account when 
determining medication regimes. Thus, for example, the researchers will often determine the 
maximum dose that can be given to a patient without risk, even though in many cases a lower dose 
provided over a longer time period may be of more benefit. 

35 Finally, it is difficult to predict the effect of the medication in the long term to anything greater than 
a very general degree, thereby further complicating the possibility of determining medication 
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Tegimes. ^ 

All of these issues are further compHcated by the profound non-linearity and complexity in ttie 
dynamics of many diseases, including disease progression and adaptation over time. 

5 

Summary of the Present Invention 

In a first broad form the present invention provides a method of determining a treatment program 
for a subject, the method including: 

a) Obtaining subject data, the subject data representing the condition; 
10 b) Using the subject data and a model of the condition to determine system values representing 

the condition; 

c) Determining one or more trajectories representing the progression of the condition in 
accordance with the model and the determined system values; and, 

d) Determining a treatment program in accordance with the determined trajectories, 

15 

The subject data typically represents the a medical condition for the respective individual, in which 
case the method includes determining trajectories representing the progression of the medical 
condition within the individual. 

20 Each system value typically represents a quantity obtained for ttie measurement of a respective 
attribute of the condition, the system values including: 

a) State variable values representing rapidly changing attributes; and 

b) Parameter values representing slowly changing or constant attributes. 

25 The method usually includes determining control variable values, the control variables representing 
attributes of the condition that can be externally controlled. 



The model generally includes one or more model equations representing the condition, the method 
including detemmiing one or more subject equations in accordance with the model equation(s) and 
30 the system values. 



The method of determining the treatment program typically includes: 

a) Evaluating the behaviour of trajectories representing solutions of the subject equations; and, 

b) Determining one or more control programs, each control program including a sequence of 
35 control variable values that result in trajectories having desired behaviour. 



' wo 2004/027674 



-3- 



PCT/AU2003/001232 



Typically the method includes determining a set of target points, the target points including stable 
points for the subject equation(s). 

The desired behaviour usually includes at least one of: 
5 a) The trajectories are acceptable; 

b) The trajectories do not move away from the target points; and, 

c) The traj ectories finally approach the target points. 

The method generally includes determining the solution trajectories to be acceptable if they are: 
10 a) Non-chaotic; and, 

b) SufBciently smooth. 

The method of evaluating the behaviour of the trajectories can include: 

a) Determining regions of control variable. and/or parameter values for which the trajectories 
15 are chaotic; and, 

b) Determining ranges of the control variable and/or parameter values for which the 
trajectories can be made non-chaotic, or otherwise stabilised. 

In tiiis case, the method typically including determining one or more control progran^is in 
20 accordance with the determined ranges. 

The method can include using a Liapunov function to determine the one or more control programs, 
although other techniques may also be used. 

25 In the case in which a Liapunov function is used, the method preferably includes: 

a) Defining a Liapunov function for which the gradient defines trajectories moving towards 
the target points; 

b) Defining constraints on the control variable values; and, 

c) Determining control variable values that result in trajectories travelling down the gradient 
30 of the Liapunov fimction in accordance with the constraints. 

The constraints generally include limits on the treatment that can be provided to the subject. 



35 



The method typically includes determining a treatment in accordance with one or more of the 
determined control programs by: 
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a) Viewing a representation of the trajectories, the representationdncluding an indication of 
the chaotic regions; and, 

b) Selecting a control program in accordance with the represented trajectories. 

5 The method may also further include: 

a) Determining one or more Nature values, the Nature values being quantities of Nature 
parameters and/or variables representing attributes of the condition that will cause the 
condition to progress in an undesirable manner; and, 

b) Modifying the subject equations to incorporate the one or more Nature values; 

10 c) Evaluating the behaviour of modified trajectories representing solutions of the modified 

subject equations; and, 
d) Performing at least one of: 

i) Determining one or more control programs, each control program including control 
variable values that result in modified trajectories having desired behaviour; and, 
15 ii) Determining one or more undesired programs, each undesired program including 

Nature values that result in modified trajectories having undesired behavioiar. 

Typically, in this case, the Liapunov function is modified to take into account the modified 
trajectories. 

20 

The method generally includes determining a set of undesired points, the undesired behaviour 
including at least one of: 

a) The modified trajectories are unacceptable; 

b) The modified trajectories do not move away firom the undesired points; and, 
25 c) The modified trajectories finally approach the undesired points. 

The method typically includes: 

a) Defining a second Liapunov function for which the gradient defines modified trajectories 
moving towards the undesired points; 
30 b) Defining constraints on the Nature values; and, 

c) Determining Nature values that result in modified trajectories travelling down the gradient 
of the second Liapunov function in accordance with the constraints. 



35 



The method usually includes determining the treatment program in accordance with control 
programs and the Nature programs by: 

a) Determining starting points having modified trajectories for which control programs exist; 
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b) Detennining starting points having modified trajectories for which Nature programs exist; 

c) Viewing a representation including at least one of: 

i) The modified trajectories; 

ii) The starting points; and, 

iii) The chaotic regions; and, 

d) Selecting a control program in accordance with one or more represented trajectories. 

The method usually includes detennining parameter values by: 

a) Detennining a partial set of system values firom the subject data; 

b) Selecting one or more models, each model including a one or more equations representing 
the effect of a condition on an individual; 

c) Attempting to determine a complete set of system values in accordance with the detennined 
partial set of system values and the respective equations; and, 

d) Selecting a model in accordance with the detennined complete set of system values. 

Hie method of attempting to detennine the complete values typically includes: 

a) Determining a candidate set of system values in accordance with the detennined partial set 
of system values and the equations; and, 

b) Comparing the candidate set of system values to at least one of: 

i) The partial set of system values; and, 

ii) Predetermined thresholds; and, 

c) Selecting the model in accordance with the result of the comparison. 

The subject is typically a patient, and may be human or non-human. Alternatively, the subject may 
comprise in vitro samples, such as cells, or the like. 

The treatment is typically the administration of medication, but may altematively, may be the 
detennined modification of any other factor effecting the condition, such as the requirement for the 
subject to receive predetermined nutrition, exercise, or the like. 

In a second broad fonn the present invention provides apparatus for detennining a treatment 
program for a subject, the apparatus including a processing system adapted to: 

a) Obtaining subject data, the subject data representing the condition; 

b) Using the subject data and a model of the condition to determine system values representing 
the condition; 



wo 2004/027674 



-6- 



PCT/AU 2003/00 1232 



c) Determining one or more trajectories representing the progression of the condition in 
accordance with the model and the determined system values; and, 

d) Determining a treatment program in accordance with the determined trajectories. 

In general, the processing system is adapted to perform the method of the first broad form of the 
invention. 

In a third hroad form the present invention provides a computer program product for determining a 
treatment program for a subject, the computer program product including computer executable code 
which when executed on a suitable processing system causes the processing system to perform the 
method of the first broad form of the invention. 

In a fourth broad form the present invention provides a method of determining system values 
representing a subject's condition, the mettiod including: 

a) Obtaining subject data, the subject data representing the condition; 

b) Determining a partial set of system values from the subject data, each system value 
reinresenting a quantity obtained by the measurement of a respective attribute of the 
condition; 

c) Selecting one or more models, each model including a one or more equations representing 
the effect of a condition on an individual; 

d) Attempting to determine a complete set of system values in accordance with the partial set 
of system values and fee respective equations, for each model; and, 

e) Selecting a model in accordance wife fee determined complete set of system values. 

The system values typically include: 

a) State variable values representing rapidly changing attributes; and 

b) Parameter values representing slowly changing or constant attributes. 

The mefeod of atten:q)ting to determine fee complete values typically includes: 

a) Determining a candidate set of ^stem values in accordance wife fee determined partial set 
of system values and fee equations; and, 

b) Comparing fee candidate set of system values to at least one of: 

i) The partial set of system values; and, 

ii) Predetermined thresholds; and, 

c) Selecting fee model in accordance wife fee result of fee comparison. 
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In a fifth broad form the present invention provides apparatus for determining system values 
representing a subject's condition, the apparatus including a processing system adapted to: 

a) Obtaining subject data, the subject data representing the condition; 

b) Determining a partial set of s: .:em values from the subject data, each system value 
5 representing a quantity obtained Vy the measurement of a respective attribute of the 

condition; 

c) Selecting one or more models, each model including a one or more equations representing 
the effect of a condition on an individual; 

d) Attempting to determine a complete set of system values in accordance with the partial set 
10 of system values and the respective equations, for each model; and, 

e) Selecting a model in accordance with the determined complete set of system values. 

Typically the apparatus is adapted to perform the method of the fourth broad form of the invention. 

15 In a sixth broad form the present invention provides a computer program product for determining 
system values representing a subject's condition, the computer program product including computer 
executable code which when executed on a suitable processing system causes the processing system 
to perform the method of the fourth broad form of the invention. 

In a seventh broad from the present invention provides a method of determining the effectiveness of 
20 treatment provided to a subject, the method including: 

a) Obtaining subject data, the subject data representing the condition; 

b) Using the subject data and a model of die condition to determine system values representing 
the effect of the condition; 

c) Providing treatment to the subject; 

25 d) Repeating steps (a) and (b) to determine modified system values; 

e) Comparing the parameter values and the modified system values; and, 

f) Determining the effect of the treatment in accordance with the results of the comparison. 

The method of determining the system values is preferably a method according to the fourth broad 
30 form of the invention. 

In an eighth broad form the present invention provides apparatus for determining tiie effectiveness 
of treatment provided to a subject, the apparatus including a processing system adapted to: 
a) Obtaining subject data, the subject data representing the condition; 
35 b) Using the subject data and a model of the condition to determine system values representing 

the effect of the condition; 
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c) Providing treatment to the subject; 

d) Repeating steps (a) and (b) to determine modified system values; 

e) Comparing the parameter values and the modified system values; and, 

f) Determining the effect of the treatment in accordance with the results of the comparison. 

5 

Typically the apparatus performs the method of the seventh broad form of the invention. 

In a ninth broad form the present invention provides a computer program product for determining 
the effectiveness of treatment provided to a subject, the computer program product including 
10 computer executable code which when executed on a suitable processing system causes the 
processing system to perform the method of the seventh broad form of the invention. 

Brief Description of the Drawings 

An example of the present invention will now be described with reference to the accompanying 
1 5 drawings, in which: - 

Figure 1 is a flow chart outlining the process of the present invention; 

Figure 2 is a schematic diagram of an example of a system for implementing the present invention; 
Figures 3 A and 3B are a flow chart detailing the process of determining a medication program; 
20 Figures 4A and 4B are a flow chart detailing the process of determining the solution trajectories for 
the method of Figures 3 A and 3B; 

Figure S is an example of a representation generated during the process of Figure 4B; 
Figure 6 is a flow chart detailing the process of determining a medication program including 
accounting for Nature values; 
25 Figures 7 A and 7B are a flow chart detailing the process of determining the solution trajectories for 
the method of Figure 6; 

Figure 8 is an example of a representation generated during the process of Figure 7B; 

Figure 9 is a schematic diagram is an alternative example of a system for inq^lementing the 

invention; 

30 Figure 10 is a schematic diagram of one of the end stations of Figure 9; 

Figure 1 1 is an example of a convergence map for a range of values of P against R, for the prophetic 
example outlined in Appendix A; 

Figure 12 is an example of a second set of convergence maps for a range of values of x(0) against P 
for fixed values of R, for the prophetic example outlined in Appendix A; 
35 Figure 13 is an example of a number of levels of magnification for a selected one of the 
convergence maps of Figure 12; 
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Figure 14 is an example of a set of target points for the prophetic example outlined in Appendix A, 
shown in greyscale and false colour, and. 

Figure ISA is an example of candidate solution trajectories for the prophetic example outlined in 
Appendix A; 

5 Figure 15B is an example c: alternative solution trajectories for the prophetic example outlined in 
Appendix A; and. 

Figure 16 is an example of candidate solution trajectories for the prophetic example outlined in 
Appendix A using a double Liapimov function. 

1 0 Detafled Description of the Preferred Embodiments 

An overview of the methodology utilised by the present invention to determine a program for a 
subject having a condition, or to detOTiiine parameter values representative of the subject's 
condition will now be described in outline with reference to Figure 1. 

15 In particular, whilst the techniques are applicable to any subject, or any condition, the techniques 
are ideally suited to determining medication programmes for patients having medical conditions, 
which will now be described in more detail. However, it will be appreciated that the techniques 
may be extended to any subject, and any condition, 

20 In any event, as shown at step 100 the first stage is for subject data to be obtained. 

Typically, the form of the subject data will vaiy depending on the nature of the subject's condition. 
Thus, for example, if it is suspected that the subject is suffering from a medical condition such as 
Parkinson's Disease, the Subject data will typically include information such as MRI brain scans, 

25 indications of brain dopamine levels, or tfie like. However, it will be appreciated that in other 
conditions, different subject data will be collected. Thus for example, in the case of a condition 
such as leukaemia, the subject data will typically include an indication of red and white blood cell 
counts, or the like. Alternatively, the subject may be an athlete, in which case, the condition may be 
the relative level of fitness for a race, with the subject data relating to indicators of the physical 

30 condition of the athlete. 

In any event, the subject data is xised at 110 to determine equations, variable values, and/or 
parameter values to represent the condition. In particular, the parameter and/or variable values can 
be inserted into equations representing the condition to determine one or more equations 
35 representing the effects of the condition on the respective subject 
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This process can be performed analytically or via numerical analysis of the subject data. 
Furthermore, aspects of this can be performed manually, although typically as the calculations are 
complex, a computational aid, such as processing system is used, as will be described in more detail 
below. 

In any event, once the parameter values, variable values and/or equations have been determined at 
step 1 10, it is then possible for the behaviour of one or more solution trajectories to be determined 
in accordance with the equations, variable values and/or the parameter values at step 120. In 
particular, with fee eq uations effectively modelling the condition of the subject, the solution 
generated by solving the equations will represent potential routes of progression of the 
c ondition within the individual. Thus, each solution tojectoiy will effectively model how the 
condition will potentially develop. ~" 

At step 130, the equations and/or the solution tiajectories are jised to -deterniine a medica^^^ or 
1 5 trt atiaaent regim e. Thusfor example, by considering different tirajectories that can be followed from 
*® ^"y ®5*!!.^SESfll^!!?^'J!i^^^ that may result in an inq)fovement, or at 

leasta.redtttai«nuaJlisjra*e^ of the subject's condition. ' 

( The. selected solution trajectories are then used to determine a medication or treatment regime. 
20 ^ Thus, for example, this may be achieved by determining the level and type of medication that may 
)he provided to the subject to allow one of the selected tiajectories which is more beneficial than the 
/ others to be followed. Alternatively, this may be used to modify the subjects lifestyle, such as diet 
I and exercise regime, to thereby improve the subject's condition, again in accordance with an 
V^propriate trajectory. 

25 

At step 140 it is possible to provide medication to the subject in accordance with the determined 
regime. 



The method can then return to step 100 allowing the process to be repeated. It wiU be appreciated 
30 that this can be performed for a number of reasons. 

Firstly, the determined parameter values, variable values, and/or equations may only be accurate 
over a short duration of time. Thus, as the condition progresses, it may be determined that the 
condition is no longer following the selected h^jectory, for example, if the determined solution 
35 trajectory is inaccurate. This may occur for example, because progression of the condition causes 
an alteration in flie model equations such that the solution trajectories determined above are only 
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accurate for the current subject condition. Accordingly, as the condition progresses, new equations, 
variable values, and/or parameter values, and hence new trajectories may need to be calculated to 
reflect the new subject condition. 

5 Secondly, new m^ chods of treatment may become available making more advantageous solution 
trajectories accessible. 

Thirdly, it may be determined tiiat the predicted solution trajectories are not being followed, for 
example indicating that the actual condition of the subject is different to that originally determined. 

10 

Fovirthly, it may be determined that the original parameter values determined are inaccurate, for 
example due to errors in the original subject data. 

It will be appreciated that other reasons may also exist 

15 

In any event, as the patent's condition progresses, either improving or worsening, it may be 
necessary to recalculate the parameter values and/or equations to allow revised medication regimes 
to be produced. 

20 This may also be performed to re-evaluate the parameter values to determine if the selected solution 
trajectories are being followed or to check if the subject's condition has improved or worsened. 

The present invention may also be used in an alternative manner as shown at steps 150 to 170. In 
this example, the subject's equations, variable values, and/or parameter values are calculated at step 
25 110. Medication and/or treatment is then proposed for and provided to the subject at step 150. 

At 160, the changes to the parameter values, variable values, and/or the equations that the provided 
medication has caused are determined. This allows the effects of the medication and/or treatment to 
be assessed at step 170. 

30 

Accordingly, this can allow drug companies or the hke to physically quantify the success of various 
forms or medication by determining actual parameter values representing the change in the subject's 
condition after being treated. This helps overcome the usual subjective assessment in which a 
subject often has to simply indicate whether they feel better after being provided a drug which does 
35 not necessary provide the drag company with useftil information regarding whether the dmg has 
caused an improvement in the individuals medical condition. 
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Similarly, in the case of performance animals, this can be used to determine the effect of training 
and diet to improve their race fitness and health. This is especially useful when dealing with 
r.nimals in which it is even harder to obtain feedback regarding the animals current state of well 
being. 

A person skiUed in the art will appreciate that aspects of the above outlined procedure may be 
performed manually. However, in order to achieve this it will be necessary for an individual to 
perform significantly complicated mathematics in order to analyse the subject data, obtain the 
parameters and/or equations and then determine the medication and/or treatment regime, hi 
particular, some of the calculations generally can only be solved using numerical methods that are 
beyond the scope of an individual. 

Accordingly, the invention is typically performed using a processing system that is adapted to 
receive subject data, and use this to determine either a medication or treatment program, or 
parameter values that are indicative of the status of a condition. An example of a suitable 
processing Systran is shown in Figure 2. 

For the purposes of simplicity, the remainder of the description will focus on the determination of a 
medication program for a patient, but it wUl be appreciated that the techniques will equally apply to 
other forms of treatment and other conditions. 

In particular, the processing system 10 generally includes at least a processor 20, a memory 21, and 
an input device 22, such as a keyboard, an output device 23, such as a display, coupled together via 
a bus 24 as shown. An optional external interface 25 may also be provided, as will be explained in 
more detail below. 

The subject data is provided to the processing system either via the input device 22 or the external 
interface 25, and the manner in which this is achieved will depend on the nature of the subject data. 
Thus for example, in the case of the subject data being an MRI scan, the scan may be suppHed 
directly to the processing system, which is then adapted to analyse the scan to extract the required 
information. Alternatively, a medical practitioner may be required to evaluate the scan to detemiine 
information such as the total brain cell mass therefrom, with this information then being submitted 
to the processing system. 



any event, once the subject data has been received, the processing system 10 is adapted to 



• wo 2004/027674 PCT/AU2003/001232 

-13- 

execute appropriate applications software stored in the memory 21, to allow the subject data to be 
analysed to thereby determine the parameter values, variable values, and/or medication program. 

Accordingly, it will be appreciated that the processing system may be any form of processing 
5 s: 7.tem suitably programmed to perform the analysis, as will be described in more detail below. 
The processing system may therefore be a suitably programmed computer, laptop, palm computer, 
or the like. Altematively, specialised hardware or the Uke may be used. 

In any event, an example of the operation of the processing system 10 in determining a medication 
10 regime will now be described in more detail with reference to Figures 3 to 6. 

In particular, the method of determining a medication regime can be performed at two levels 
depending on the medical condition and the level of accuracy required. 

15 In the first scenario outlined in Figures 3 and 4, it is assumed at first instance that the medical 
condition is passive. Thus effectively, this assumption indicates that the medical condition will not 
respond in an unexpected mamier to the provision of medication. Thus, it is assumed that if the 
medication will counteract the effects of the condition, the condition will not then respond to 
counteract the effects of the medication. 

20 

It will be appreciated that this technique is not necessary suitable for all medical conditions but will 
in any event provide a first level of approximation that can be used in determining a medication 
regime. 

25 In any event, as shown in Figure 3 A the process begins at step 200 when subject data is obtained by 
observation of the subject. The subject data is then used to determine a partial set of state variable 
values and parameter values at step 210. The state variable and parameter values z, X 
represent quantities obtained for the measurement of respective attributes of the medical condition, 
with the partial set , being tiiose that can be determined fi-om the subject data. 

30 

Jn particular, the state variable values z are values of measurements that are subject to rapid 
changes, typically on a day to day scale. Thus, for example, in tiie case of Parkinson's Disease a 
state value may reflect the current level of dopan[une within the medicated brain, which generally is 
subject to significant fluctuations on an hoiir to hour basis. 

35 
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In contrast, the parameter values X are values of measurements that tend to remain relatively *-T 
constant on a day to day basis, having only long-term variations. These are therefore slowly 
changing values that typically will change gradually over a number of months or years. Thus, for 
example, in the ca : : of Parkinsons Disease, a parameter value may include the number of living 
5 dopaminergic brain ce' Is in the individual. It will be appreciated this number of brain cells will 
typically decay gradually over a few years and is therefore not subject to the rapid change. 

It will be appreciated that the partial set of state variable and parameter values Z^^X^ may therefore 
be determined in a number of manners. Thus, for example, the subject data may simply be provided 
10 to the processing system 10, which then extracts the partial set of state variable and parameter 

values. Alternatively however the partial set of state variable and parameter values can be ( 
detemiined manually by a user and then input into the processing system 10. 

Furthermore, a database of predetermined state variable and parameter values may be maintained 
IS for different medical conditions. The predetermined values may be based on data collected from a 
number of subjects, and therefore forms generic data regarding the progression of the medical 
condition in a sample population rather than specific data concerning the progression of the medical 
condition in a respective individual. 

20 Accordingly, this allows the medical practitioner to select a partial set of state variable and 
parameter values that correspond to the current level of progression of the medical condition in the 
individual, which is in turn determined from the subject data. In this case, as the determined values 
are representative of the medical condition in a general population, any detemiined medical ( 
program will not be specific to the individual, but will be generic to the medical condition. 

25 

However, this may provide a sufficiently accurate medication program for many cases. Otherwise, 
the partial set of state variable and parameter values are determined such that they are specific to the 
subject, allowing a subject specific medication program to be developed. 

30 In any event, having determined the partial set of state variable and parameter values z", X^ at step 
210, the processing system 10 operates to select one or more predetermined models at step 220, 



35 



In general, the processing system 10 will have a number of predetermined models stored in the 
memory 21. Each of these models will include one or more respective partial and/or ordinary 
differential equations that represent the progression of a medical condition for a subject. 
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It will be appreciated that a respective model may be provided for each condition. In addition to 
this, several different models may be provided for a respective condition. This may be necessary 
for example if the effect of a condition is different in different individuals, or if the effect of a 
condition varies depending on the current state of the condition within the subject. Thus, for 
5 example, the effect of Parkinson's Disease on an individual will tend to be very different at an early 
stage, compared to at a later stage when the condition has progressed further. Accordingly, 
conditions such as Parkinson's Disease will tend to have a significantly different set of ordinary 
dififerential equations during initial stages of the condition when compared to the situation in which 
the condition has progressed significantly within the individual. 

10 

Accordingly, the processing system 10 selects one or more of the predetermined models, which on 
the face of it may model the condition of the subject. Thus for example if the determined condition 
is Parkinson's Disease, then the processing system 10 can simply select one or more different 
Parkinson's Disease models. 

15 

In order to make the selection, the processing 10 can perform the selection on the basis of input 
commands provided by a user. Thus for example a trained medical user may provide their own 
assessment of the condition. Alternatively however the processing system 10 can be adapted to 
select models in accordance with the partial set of state variable and parameter values z^,X^, 
20 provided at step 210. 

Li any event in steps 230 and 240 the processing system 10 attempts to fit the determined partial set 
of state variable and parameter values into one of the models. 

25 In order to achieve ttiis at step 230 the processing system 10 atteix^>ts to construct a complete set of 
state variable values z and parameter values X from the partial set of state variable and parametQ:s 
values z'*,^'' and the respective model equations. Thus, the equations will typically specify various 
relationships between the different state and parameter values, and accordingly, the processing 
system uses these relationships to determine values for the state variable and parameter values z, X 

30 that have not been determined from the subject data. 



At step 240, the processing system 10 evaluates each model in accordance with the determined 
complete set of state variable and parameter values z, X. 
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Thus for example, the processing system 10 will compare the equations for each model to determine 
whether the determined state variable and parameter values z, A. are self consistent, and to determine 
if these fit with the partial set of state variable and parameter values z^^X^ determined from the 
subject data. 

5 

If it is determined that none of the models are acceptable at step 250 the process returns to step 220 
to allow a different model to be selected. Steps 220 to 250 are then repeated until a suitable model 
is selected. 

10 It will be appreciated that if no suitable model is available, then it may be necessary to determine a 
new model. This is achieved by performing a numerical analysis of state variable and parameter 
values for a subject suffering from the medical condition. In particular, variations in the parameter 
and state variable values over time are determined, and then used to generate partial or ordinary 
differential equations that best model the changes in the values over time. It will be realised that the 

15 developed equations can then be used to predict further changes in parameter values, thereby 
allowing the differences between predicted and obtained values to be used in improving the model. 

In any event, once a suitable model has been accepted at step 250, the process moves on to step 260 
to allow the processing system 10 to determine a subject specific model, of the form: 
20 dz/dt=f(z,uAt) (1) 

where: 

z is a state vector formed from the state variable values such that z € A c SR" 
A is a set of all possible state variable values 
25 u is a control vector formed from control variable values such that u € U ^^R^ 

U is a set of all possible control variable values 

X is a parameter vector formed from the parameter values such that \e A c SI** 
A is a set of all possible parameter values 
t is time 

30 

The subject's specific model is obtained by including the state variable and parameter values into 
the model of equations determined at step 220 to 250 above. In addition to this, the equations will 
also include a time dependency represented by time t, and a control aspect defined using the control 
vector u. 
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The control vector u represents various external factors that can be used to influence the progression 
of the condition, such as the application of medication, or the like. At this stage, the control vector 
is generalised and does not include specific values. 

5 Once the equations have been determined, this allows the system equations to be used to generate 
solution traj ectories <f>, such that 0 (z, u, \ t) c 91 ° . . 

At step 270, the processing system 10 operates to locate reachable equilibrium or other desired 
points for the subject specific model. The equilibrium points are given by points at which the rate 
10 of change of the subject's specific model with respect to time is zero, such that: 

dz/dt=0 (2) 

This is achieved by contiputing an initial state vector for which e A , and for which 
15 au*^ eU, A-'^ € A suchthat: 

lim^^T^z***^*^ ^^^'O-^z'^e^^fOTsome T<oo (3) 

where: 

f( 2 * , u ^ A ^ , t) s 0, f or all t > T. 

20 

Thus, the processing system 10 operates to calculate initial state variable values z° for which 
solution trajectories <p can be steered towards the equilibrium or other desired points as the time t 
tends towards a sufficiently large value. These solution trajectories <p represent progressions of the 
condition that lead to equilibrium or other desired points. 

25 

It will be appreciated however that the trajectories q) may not be trajectories that can be followed in 
real terms by the subject as they may require excessive levels of drug dosage, or the like. In 
addition to this, not all of the equilibrium points would represent desirable outcomes for the subject 
as some of these equilibrium points may correspond to situations in which the condition is 
30 significantly worse tfian at present. Thus, for example, a equilibrivmi point will typically occur 
when the subject is dead, which is obviously an undesirable outcome. 



Accordingly, at step 280 the next stage is for the processing system to determine those solution 
trajectories that could potentially be followed by the subjects in real life. In order to do ttiis, the 
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processing system 10 contains a constrained range of control variable values and parameter 
values for which solution trajectories are acceptable from a particular initial position. 
Accordingly, this corresponds to determining trajectories emanating from an initial state for 
which the trajectories can be made acceptable. 

5 

In this sense, acceptable solution trajectories are solution trajectories that are non-chaotic and which 
are sufficiently smooth so that they do not adversely affect the subject. In particular, it is generally 
not advisable for a subject's condition to be induced to follow a rapidly changing trajectory as this 
will result in rapid changes in the subject's condition which will generally be detrimental to the 
10 subject's overall health. 

.( 

In order to achieve this, the processing system operates to test for stability as the parameter values 
'k vary. 

15 La particular, the processing system 10 operates to determine parameter values X for which chaos 
emerges on a subset of all possible state and paramet^ values, given by: 

A*®A* c:A®A (4) 

where: 

20 0 is \ised to denote the Cartesian product of subspaces 

Once a position has been determined for which chaos occurs, the processing system 10 operates to 
compute eigenvalues ^ (z, u, \ t) and constrained control and parameter values u , such that the ^ 
eigenvalues ^(z, u*^ , , t) stabilise the system. This is also used to designate avoidance sets: 

25 A=A^®A^ (5) 

where: 

A are the avoidance sets 
A^^®A'='> cA. 



30 Thus, this process allows the extremities of chaotic regions in state space to be determined. This in 
turn allows the chaotic state space regions to be avoided and/or responded to by the control 
programs that are subsequently determined, as will be described in more detail below. 
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In particular, the processing system 10 determines sets U*" , A*" of constrained control and parameter 
values u%X% and uses these sets to determine a region of controllability R*" c: A, in state space. 
This is performed such that for z^^eR", the eigenvalues ^(z^ u^ l^ t) stabilise the system trajectories 
cp into non-chaotic behaviour. 

5 

Accordingly, this allows the processing system 10 to determine ranges for the control variable 
values u*^ and parameter values X"" for which the solution trajectories 9 are non chaotic within the 
state space. 

10 Having detennined regions of chaos at step 290, and then used this to determine state space regions 
where acceptable solution trajectories cp exist, the processing system 10 operates to determine a 
family of control programs p*^ for which acceptable solution trajectories <p will be obtained. In 
particular, the control programs p** are based on the ranges of the constrained control variable 
values u** for which acceptable solution trajectories exist 

15 

Thus, given any initial state vector 6 R° , the control programs p*^ can be determined such that 
the solution trajectories 9 are non-chaotic: 

p'^(2A%t):R'->U" (6) 

where: 

20 4>iz^,p'' (zA%t),X%t) is non-chaotic. 

Having determined a family of control programs p*' this allows a specific control program to be 
used to determine a medication regime at step 310. This can be performed in a number of different 
ways. 

25 

Thus, for example, the processing system 10 could be adapted to generate all possible solution 
trajectories 9 that are acceptable and then present these to a user allowing a user to select a solution 
trajectory 9 which appears appropriate. However, this will typically be computationally very 
expensive and time consuming and therefore would be undesirable particularly if the system is used 
30 on a day to day basis. 



Accordingly, Figures 4A to 4B show a particular technique for achieving this using a Liapunov 
function. 
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As shown at step 400 in Figure 4A, this is achieved by having the processing system 10 select a set 
Xi of desirable target points for the subject's specific model dz/dt = f(z, u, \ t). 

5 In general, the target points will typically correspond to selected points in he neighbourhood of the 
equilibrium points determined at step 270 above. Alternatively, points that are otherwise 
determined to form stable sets may also be used. However, it will be appreciated that this is not the 
only criteria which would be used to select target points z^ and in particular, any equilibrium or 
stable points corresponding to a deterioration when con^ared to the subject's current condition will 
10 generally be excluded. However, this may not be the case if the subject has a temiinal illness for 

nample. In this case, whilst the stable points may be worse than the subject's current condition they ( 
may be the best that the subject can hope to achieve. 

In any even^ the selection of the target points z^ may be achieved auton:mtically by the processing 
15 system, in accordance with predefined criteria. Alternatively, the target points may be selected 
by a medical practitioner, for example, by evaluating the suitability of the stable points as target 
points 2^ . 

At step 410 the processing system determines a Liapunov function V^. The Liapunov function is 
20 based on desired qualitative behaviour to be imposed on system trajectories and it will therefore be 
appreciated that a number of different Liapunov functions may be selected for any on particular 
scenario. ^ 

It will thorefore be something of an art for a user of the system to select tfie most advantageous 
25 Liapunov function, the first time a respective set of equations are solved. However, in general even 
though the specific equations of the medical condition will vary between subjects it is generally 
possible for similar equations to use similar Liapunov functions when solution trajectories are being 
determined. Accordingly, any Liapunov functions used will be stored in the memory 21 so that a 
database of suitable Liapunov functions is created. Subsequently, Liapunov functions can be 
30 selected firom the database as required. 

In any event, the Liapunov function is determined such that: 

V, :A®A->5R, (7) 



35 



In any event, at step 410 during the creation of the Liapunov function the processing system will 
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operated to define constraints on the parameter values X and the state variable values z. These are 
determined based on the regions of chaos determined at step 280 above. 

The safety constraints are determined such that X^A\ where A is a bounded set such that 

5 S A . Thus, this typically involves using a sub-set of the set of constrained parameter values 

A*^, determined in step 280 above, as will be outlined in more detail below. 

Furthermore, a closed set of target points is determined such fliat: 

= {z| V,(z,X) < C,}VX € A\ (8) 

where: 

C, > 0, such that 2^ e . 

At step 420 the processing system 10 operates to determine constraints on the control variable 
values u. The constraints are generally defined in accordance with (he constrained control variable 
values u^ determined in step 280 above, although additional constraints may also be applied. 

Thus additional constraints u^, u"" may also be defined based on limits on the control of the 
condition that can be provided to the subject As control is usually provided in the form of 
medication, it wQl be appreciated that the constraints are generally focused on limits for flie 
provision of the medication to the individual. This may be based on factors such as maximum drug 
dosages, maximum drag dosage rates, maximum drag absorption rates by the body, costs, or the 
like. 

However, in other conditions, other factors n^iy also be influential. This may include for exanqjle 
obtaining regular exercise in the case of high cholesterol levels or heart disease. In this case, 
restraints will be defined in accordance with exercise regimes proposed for the individual. From 
this it will be appreciated that whilst medication generally refers to the provision of drags or other 
pharmaceutical products to the individual, medication will also generally cover any form of 
treatments of the condition including for example, dietary requirements, exercise, or the like. 

In any event, the constraints are defined such that u e [u'jU*^] ^U. 



At step 430 the processing system 10 then determines one or more control programs p*, which 
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effectively control variable values u that are constrained to provide a desired solution trajectories <p. 
The control programs p* are generated to desired solution trajectories that: 

• Are acceptable; 

• Do not move away from tlie target point set ; and 

• Finally enter the target point set x, . 

Thus the at this stage, the processing system 10 efifectively operates to exclude solution trajectories 
that either move away from the target point set Xj , do not finally enter the target set Xj , or are not 
acceptable. 

As described in some detail above, acceptable solution trajectories are trajectories that are non- 
chaotic and are sufficiently smooth. Accordingly, it will be appreciated that ttie desired solution 
trajectories vidll effectively be a sub-set of diose trajectories determined to be acceptable at step 280 
above. 

Thus, the constrained parameter values X* must at least be a sub-set of the constrained parameter 
values X^^/l with the control program p* also being a sub-set of the control variable values, u^elT. 
This. is required to ensure that the desired solution trajectories are acceptable, as well as ensuring 
that they do not move away from the target point set and finally enter the target point set x, . 

Expressed mathematically, the processing system operates to determine the control program p*: 

p*(A-*):A\x,-^[u-,u-h], (9) 

where: 

A\xi denotes the domain of state vectors that fall outside the target set Xi 
X* are the constrained parameter values determined such that A.* e A^ 

and, for all t>0, p*(X*) attempts to achieve either: 

dVj /dt< 0 (inqjrovement as the trajectories move closer to the target point set x^ with 
respect to the contours of Vj (z) ); or at least, 

dVj /dt=0 (no worsening as flie trajectories do not move away from target point set Xj 
with respect to the contours of V, (z) ). 

In any event, once one or more control programs which result in acceptable solution trajectories are 
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determined, the processing system 10 generates candidate solution trajectories cp for each control 
program p*. This is performed by the processing system 10 at step 440. 

In particular, the solution trajectories are determined in accordance with the control program p* and 
5 the constrained parameter values X*. In addition to this, the solution trajectories are determined in 
accordance with the current state of the subject, which defines an initial state vector z°, such that the 
solution trajectories are as follows: 

<Kz',P*.^*.t) (10) 

10 At step 450 the behaviour of each candidate solution trajectory 9 is assessed. This may be achieved 
in a variety of manners depending on tiie implementation. 

Thus for exaiiq>le, the processing system 10 can be adapted to compare the solution trajectories to 
predefined criteria. This may include for example conqparing end points of each candidate solution 
IS trajectory <p to see which end point is closest to a selected one of the target points z^, or to 
determine which solution trajectory end point would be most beneficial for the subject. Li addition 
to this, the processing system may be adapted to determine which candidate solution trajectories 
<p remain furthest fi:om the chaotic regions at all time, or the like. 

20 Alternatively, the candidate solution trajectories may be assessed by a user of the system. 

In order to achieve this, the processing system 10 can be adapted to generate a representation on the 
display 23. An example representation is shown in Figure 5. 

25 As shown, the representation would in general be a two dimensional graphical representation of the 
state space showing the solution trajectories q>, and regions where solutions to the subject specific 
equations are chaotic C. The representation could also include examples of the target point set Xj , 
and optionally the initial state variable values z°. 

30 This will allow a medical expert to select a solution trajectory fliat would appear to be most 
beneficial to the subject. Thus, for example, the trajectory cpi leads from the initial start vector to 
the desired target point set, and would tiierefore be preferable to other solution trajectories 92^ 93 
which do not lead to the desired target point set, or which do not start from the initial start vector 2°. 
However, it will be appreciated that combinations of solution trajectories may also be selected as 

35 shown at 94, 95. 
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in addition to this it is also possible to select a desired trajectory cp, or modify the control program 
p* to determine alternative trajectories cp based on either the qualitative behaviour of the trajectories 
0(z^, p*, X*, t) or through quantitative optimisation. This can be achieved for example via 
5 analytical methods, such as Euler equations or via adaptive numerical methods, such as through the 
use of genetic algorithms. 

At steps 460 it is determined if a desired solution trajectory (p has been determined. If so, the 
processing system 10 operates to provide an indication of the desired solution trajectory (p and tiie 
10 associated control program p* at step 470. 

This will include a numerically computed solution trajectory <^z^, p*, X*, t) and final control 
program p*, as well as a region of finite controUabiUty c: A determined such that 
Vz*' € i?^,lim,^3. <Kz^ p*. A.*, t) ^ 2' € Ti , for some T>0. 

15 

Otherwise the processing system modifies one or more of the Liapunov function V, , the constraints 
u , and the target point set before r^eating steps 400 onwards as required at step 490. This 
is represented by the link back to 4A as shown. 

0 Again, modification of the constraints and the target points set can be achieved in a number of 
ways, such as via analytical or adaptive numerical methods. Alternatively, this can be achieved 
manually via consideration of the constraints and the target points set. 

Thus, for example, a medical practitioner can consider the constraints determined based on the 
> treatment, such as the medication available, and consider whether this can be adapted. Thus for 
example, can alternative drugs, or higher doses be used. Similarly, if the desired target points t, 
cannot be reached, can other slightly less desirable target points be reached. 

As described above, this example relies on the condition being passive. However, in many case, the 
) condition is not passive and instead responds to medication or the Uke, to attempt to counteract the 
effects of the medication. Thus, for example, if medication is provided to a subject, the subjects 
condition may fight against the action of the medication by altering it's behaviour patterns in an 
attempt to counteract the effects the medication. Thus, for example, in the case of immune 
deficiency conditions, the provision of drugs to kill a virus may result in an initial increase in viral 
levels as the virus reproduces in an attempt to counteract the effect of the drugs. 
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Alternatively, steps 260-460 may have revealed that the system behaviour depends on some state 
variable or parameter value, which is measurable only up to a margin or certainty. The worst case 
scenario inherent in this margin of uncertainty must be assessed now. 

5 

In order to account for this, a modified method can be used which is classed as a "Game against 
Nature" type problem. 

Again, this form of problem can be dealt with in a nvunber of ways. However, a specific example 
10 will now be described with reference to Figures 6 to 8 in which Nature values w are used to 
represent the action of the condition to counteract the effects of any medication. This may be 
achieved by using a Nature vector w formed ftom Nature parameter values such that w e S cz SR" . 

In particular, •'Nature" is taken to be a sentient opponent that calculates the best possible choice of 
15 values for opposing the efforts of the clinician witirin the determined constraints. Accordingly, the 
Nature vector w will be designed to take into account factors such as system noise imcertainty, 
instrument uncertainty, dmg admiidstration noise and condition progression as well as other 
possible sources of uncertainty or pathology. 

20 The Nature values are therefore quantities of Nature parameters and/or Nature variables 
representing attributes of the medical condition that will cause the medical condition to progress in 
an undesirable maimer. Thus, it will be appreciated tiiat the Nature parameters are slowly changing 
attributes, with the Nature variables being rapidly changing attributes. 

25 Accordingly, at step 500, the processing system detennines various Nature parameter and/or 
variable values. 

The specific Nature values may be determined by a medical user, or the like, in accordance with 
observations of the subject and/or condition. Alternatively, Nature values may be selected from 
predetermined Nature values determined for specific conditions. 

30 

At step 510, the processing system 10 utilises the Nature values w to determine a modified subject 
specific model: 

dz/dt = f(z,u.w,X,t) (11) 
35 Again these system equations generate solution trajectories <ji(z°, u, w, t) c SR" , 
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As in the previous example, the processing system 10 operates to locate equilibria aiid qualitative 
system behaviour for the subject's specific model dz/dt, at step 520. 

The processing system 10 flien determines a range of control variable values u"""^ and parameter 
values X*^'' for which solution trajectories 9 are acceptable. 

Again, this will involve testing the solution trajectories (p to find regions of chaotic behaviour in 
state space, in a manner similar to that described above with respect to Figure 3B. The processing 
system 10 will then compute eigenvalues ^(2, u, w, X, t), and control variable and parameter 
vectors u^, X*^^ such that ^ (z, u^, w, X''^ t) stabilises flie system for Vw e 3 . The processmg 
system 10 uses these values to determine the avoidance sets A** such that A** ® A*"** c A** . 

Having determined state space regions for which acceptable solution trajectories exist despite the 
effects of the Nature vector w, at step 540, flie processing system moves on to step 550 to determine 
a family of control programs p*^ in accordance with the determined ranges of control variable 
values u*^ and parameter values X"^. 

In particular, the processing system 10 computes a region of controllability R*^ in state space, such 
that given any initial state vector z° eR*^, the control programs p*''" are based on the control 
variable values u'^^Cz^, w, X'^^ t) such that solution trajectories cp(z^ u*=^(z^ w, X"^ t) w, X*^^ t) are 
non-chaotic for Vw e H , despite the effects of the Nature vector w. 

At step 560, the femily of selected control programs p*'"' is used to determine a medication regime. 

Again, the procedures set out in Figure 6 may be performed in a number of manners. However, in 
this example, the processing system determines a second Liapunov function that is used to represent 
the effects of Nature, such as the effects of the medical condition and other noise or uncertainty 
factors. 

Thus, as shown in Figure 7A, at step 600 the processing system 10 determines a modified first 
Liapunov function Vi^ and corresponding target point set Ti^. In this case, the Liapunov function 
V,w is a modified version of Liapunov function Vj used in Figures 4A and 4B, which has been 
modified to take into account the Nature values w. 

Once this has been completed, the processing system 10 then operates to select a set T2 of undesired 
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points -f^ at step 610. These undesired points will correspond to lethal or extremely undesired 
points which maybe detrimental to the health of the subject if the subjects condition were to pass 
through these locations. 

The set of undesired points can be determined based on: 



where: 

C2>0 

This allows the processing system 10 to define a large avoidance set A = ® A^ such that 
A^'cA^cA^t^qA^cA. 

At step 620, the processing system 10 determines a second Liapunov function V2. The purpose of 
this second Liapunov function V2 is to model the operation of the Nature values w as a hostile 
entity, to thereby embody "enemy" behaviour of noise/uncertainty and condition progression. 
Accordingly, the second Liapunov function V2 determines solution trajectories representing the 
progression of flie medical condition under the influence of the Nature vector, which causes a 
worsening in the medical condition, as far as controlling or preserving the subject's health is 
concerned. These solution trajectories will hereinafter be referred to as undesirable trajectories. 

The second Liapunov function is defined such that : A <8) A — >• 9? . 

At step 630, the processing system 1.0 determines acceptable constraints on the Nature values w. 

As described above, the Nature values w are determined on the basis of factors such as the Systran 
noise, instrument imcertainty, drug administration noise, or the like. Accordingly, a nvimber of 
these values will be determined based on the manner in which subject data was obtained, the 
method of drug delivery, inaccuracies in drug administration methods and the like. 

This can be used to determine flie areas that introduce the most uncertainty or lack of confidence 
into deterniining and treating the subject's condition. Thus, for example, if it is determined that 
there is a high level of confidence associated with existing tolerances in measurements of ligand, it 




(12) 



A**eT2®A 



(13) 
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will be appreciated that massive investment into obtaining more accurate measurement tolerances 
would be relatively worthless. 

Accordingly, assessing the Nature values in itself is useful, as it will allow areas of low control 
confidence to be identified, allowing these issues to be addressed in preference to other issues 
which introduce less threat 

It any event, at step 640, the processing system 10 determines one or more Nature programs p** 
which cause the undesirable trajectories acting under the influence of the Liapunov functions V2 to 
move towards the set of undesired points T2. This is performed at step 640 to allow the processing 
system 10 to determine the worst effects of changes in the medical condition. Thus, this effectively 
represents Ae circumstances in which the medical condition deteriorates or an uncertainty presents 
itself in the worst possible fashion, as far as controlling or preserving the subject's healtii is 
concerned. 

In particular, for a given initial state vector z'^eA\(x2^A'^) the processing system operates to 
determine a Nature program p* * given by: 

p**:A\(xiUA^) -^[w", (14) 

where: 

A"^ is the avoidance set such that t^ctA^ 

and, for all t>0, p** attempts to achieve either: 

dVj /dt< 0 (deterioration as xindesirable trajectories move closer to the undesired point 
set T2 with respect to the contours of Vj (z)) ; or at least, 

dVj /dt= 0 (no improvement as undesirable trajectories do not move away from the 
undesired point set T2 with respect to the contours of V, (z)) . 

Accordingly, the family of Nature programs p** represents the set of the condition's strongest 
"choices" of control strategy to resist therapy. 

At step 650 the processing system 10 operates to determine one or more control programs p*"^ for 
which the solution trajectories are acceptable and do not move away from the target point set Xjw- 
Thus, this corresponds to finding solution trajectories representing the progression of the condition 
in which the solution trajectories are non-chaotic and smooth, and after which the progression does 
not worsen, despite the best efforts of Nature. 
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In particular, for a given initial state vector z°*€A\(tiwUA^) the processing system operates to 
determine a control program p*'^ such that: 

p* ^ :A\(x„vVJA^) ->[u , u^ (15) 

5 such that: 

liint^T9(z°, p*^ w. A,, t)->z"exiw ; for T^O; (16) 

V XeA^'" 

Vwe[w",w*]; 

10 where: 

I A'^ is the avoidance set, such that T2^^ 

and, for all t>0, p* attempts to achieve either. 

dVxw/dt<0 (improvement as trajectories move closer to the target point set Xiw with respect 
15 to the contours of Viw(z)); or at least, 

dViJdMi (no worsening as trajectories do not move away from the target point set Xiw with 
respect to the contours of Viw(z)) 

Thus, the processiag system 10 effectively determines initial state functions for which a control 
20 program can be determined that steers the solution trajectories towards the desired target point set 
xiw despite the best efforts of flie disease to achieve deterioration or stalemate. The control program 
p** therefore represents the clinicians' "strongly-winning strategy", with the family of such control 
programs p*"^ representing the set of the clinicians' strongest choices of control strategy to achieve 
the desired target point set Xiw- 

25 

Accordingly, from this it will be appreciated that the processing system 10 determines one or more 
Nature programs p** and one or more control programs p*"^ which represent the strongest control 
strategy for the condition and the strongest control strategy for the cUnician respectively. 

30 At step 660 the processing system 10 operates to determine a region Wi of initial state variable 
values z*^* for which control programs p*"^ exist. Thus effectively, the processing system 10 
operates to examine the solution trajectories and determine values for the control, parameter and 
Nature values for which control programs p** exist. 



35 



Thus, the processing system 10 operates to detmnine: 
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Wi={z^* €A\(x,^uA'')such^tp**' exists.} (17) 

At step 670 the processing system 10 operates to similarly determine regions W2 of initial state 
vectors z*^** for which Nature programs p** exist 

5 

Thus effectively, the processing 10 works out initial state vectors z°** which result in imdesirable 
trajectories that follow the Nature programs p** and state vectors z°* for which the solution 
trajectories follow the control programs p*"^. 

10 The processing system then generates candidate solution and undesirable trajectories for each 
control program p*"*", and each Nature program p** at step 680. 

The behaviour of each candidate trajectory is then assessed at step 690. 

15 Again, this may be achieved in a number of ways and may be achieved for example automatically 
by having the processing system analyse the candidate trajectories and choose the solution 
trajectories that will most benefit the subject. 

However, typically this is achieved by having a professional medical staff member examine a 
representation showing the regions Wi, W2, the chaotic regions C and the solution trajectories <p as 
shown in Figure 8. 

Thus as shown, regions of chaotic behaviour as shown at C, with the regions of undesired point set 
t2 and target point set Xjw and initial regions Wi, W2 as shown. 

In use, the medical practitioner would determine the subject's initial state variable values that form 
the state vector . The medical practitioner then selects one or more trajectories fh)m fliis point 
that will ultimately allow the subject's medical condition to move towards the target point set Ti^. 

In this example, seven different trajectories are shown as (pj to 97 respectively. 

In this example only a single one of the trajectories 93 extends from the initial starting point to 
the target point set Ti^. However, fliis trajectory 93 enters a chaotic region C that can therefore 
instantly be dismissed. 
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The trajectories cp6 and (p? ultimately end at the iindesired point set X2 and again, it would be 
undesirable to follow or even approach these trajectories. 

Accordingly, it will be necessary for the subject regime to cause a subject to progress along one of 
5 the trajectories <pi, (p2 to the region Wi and then follow one of the trajectories <p.v, 95 to the target 
point $et T]w. In this example the trajectories cpi, 95 both near regions of chaos C and accordingly, 
the trajectories (()2, 94 would be selected. 

Once an acceptable trajectory is determined at step 700 tiie processing system 10 provides an 
10 indication of the desired solution trajectory and the associated control programs p'^^. It will be 
{ appreciated that in this exaniple, two solution trajectories are required and accordingly two 

corresponding control programs p**' will be required. In this instance, the first control'program will 
cause the subject's condition to progress along Ifae trajectory 92 with the second control program 
causing the subject trajectory condition to move along the trajectory 94 thereby ending within the 
15 target point set Xjw- 

Jn the event that a suitable solution trajectozy has not been determined then the processing system 
moves on to step 720 to modify one or more of the Liapunov function V, , the Liapunov function 
Vo, the Nature parameter constraints w",w'^, the control parameter constraints u**,u'*'and the 
20 target point set Xiw 

^ The processing system 10 then repeats step 600 onwards as required. 

Architectures 

25 It will be appreciated that the above method may be achieved in a number of different manners. 
Thus, for example, a respective processing system 10 may be provided for each medical practitioner 
that is to use the system. This could be achieved by supplying respective applications software for a 
medical practitioner's computer system, or the like, for example on a transportable media, or via 
download. In this case, if additional models are required, these could be made available through 

30 program updates or the like, which again may be made available in a number of manners. 

However, alternative architectures, such as distributed architectures, or the like, may also be 
implemented. 
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An example of this is shown in Figure 9 in which the processing system 10 is coupled to a database 
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- 11, provided at a base station 1. The base station 1 is coupled to a number of end stations 3 via a 
communications network 2, such as the Internet, and/or via communications networks 4, such as 
local area networks (LANs) 4. Thus it will be appreciated that the LANs 4 may form an internal 
network at a doctor's surgery, hospital, ,or other medical institution. This allows the medical 
5 practitioners to be situated at locations remote to the central base station 1 . 

Accordingly, in use the end stations 3 must be adapted to communicate with the processing system 
10 positioned at the base station 1. It will be appreciated that this allows a number of different 
forms of end station 3 may be used. 

10 

An exanq>le of a suitable end station is shown in Figure 10. As shown the end station 3 includes a 
processor 30, a memory 3 1 and an input device 32 such as a keyboard, an output device 33 such as 
a display coupled together via a bus 34, as shown. An internal interface 35 is typically provided to 
allow the end station to be coupled to one of the communications networks 2, 4. 

15 

In use, the processor 30 is adapted to communicate with the processing system 10 provided in the 
base station 1 via the communications networks 2, 4 to allow the above described process to be 
implemented. Accordingly, it will be appreciated that if the commimications network 2 is fee 
Internet, this will typically be achieved by having the base station 1 present web pages to the 
20 medical practitioner on the end station 3. 

Accordingly, it will be appreciated that the end stations 3 may be formed from any suitable 
processing system, such as a suitably programmed PC, Litemet terminal, lap-top, hand-held PC, or 
the like, which is typically operating applications software to enable data transfer and in some cases 
25 web-browsing. 

In this case, the patient data, and any other information as well as input commands, may be supplied 
by the medical practitioner via the end station 3, before being transfeired to the processing system 
10, located at the base station 1. The processing system tiien operates as described above, 
30 transferring the results to the end station 3 for presentation to the medical practitioner. 

In this case, it will be appreciated that access to the process may be controlled using a subscription 
system or the like, which requires the payment of a fee to access a web site hosting the process. 
This may be achieved using a password system or the like, as will be appreciated by persons skilled 
35 in the art. 
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Furthermore, information may be stored in the database 11, and tins may be either the database 11 
provided at the base station 1, the database 1 1 coupled to the LAN 4, or any other suitable database. 
This can also include patient data, results, determined parameters, or the like. This allows the 
medical practitioner to maintain a record of patient's medic condition history, which can be 
5 accessed as required at subsequent times. This may be used as an accurate record of the progression 
of the medical condition within the patient, and if coupled with a record of the medication actually 
provided to the patient, this can allow the effect of the medication to be determined. 

It will be appreciated that by analysing data collected for a number of patients, this will 
10 allow more accurate models to be developed in an iterative process. Statistical analysis can 
also allow additional model to be developed, for example by analysing a range of age 
groups to create age dependent models. 

Example 

15 A prophetic example of the technique as supplied to a colony of cells in the human body is 
shown in Appendix A, below. 

Variations 

It will be appreciated from the above, that the techniques can be applied to any subject, and 
20 this includes, but is not Umited to patients of hiunan or other mammaUan, or non- 
mammalian species and includes any individual it is desired to examine or treat using the 
methods of the invention. Suitable subjects that fall within the scope of the invention 
include, but are not restricted to, primates, livestock animals (e.g., sheep, cows, horses, 
donkeys, pigs), laboratory test animals (e.g., rabbits, mice, rats, guinea pigs, hamsters), 
25 companion animals (e.g., cats, dogs) and captive wild animals (e.g., foxes, deer, dingoes). 

It will also be appreciated that the techniques can be used in vitro to exancdne the reaction of 
specific samples. Thus for example, the teclmiques can be used to monitor the reaction of 
cells to respective environmental conditions, such as combinations of nutrients or the like, 
30 and then modify the combination of nutrients, to thereby alter the cells response. 

Furthermore, it will be understood that the temis "patient" and "medical condition" where 
used do not imply that symptoms are present, or that the techniques should be restricted to 
medical or biological conditions per se. Instead the techniques can be appUed to any 



wo 2004/027674 



-34- 



PCT/AU2003/001232 



condition of the subject. Hius, for example, the techniques can be applied to performance 
subjects, such as athletes, to determine the subject's response to training. Tliis allows a 
training program to be developed that will be able to prepare the subject for performance 
events, whilst avoiding overtraining and the like. 

Thus, it will be appreciated that the condition of the subject may be the cun-ent physical condition, 
and particularly the readiness for race fitness, with the treatment program being a revised training 
program specifically directed to the athletes needs. 

Thus, it will be appreciated that the techniques outlined above can provide: 

(1) a method and apparatus for determining a medication program for a patient. 

(2) a method and apparatus for determining a veterinary medication program for an 
animal suffering from a medical condition, and to a method and ^paratus for 
determining livestock parameter values representing the effect of a medical 
condition on the subject animal; 

(3) a method and apparatus for determining a non-linear control program for 
pharmaceutical, pharmacological or physiological processes, and to a method and 
apparatus for determining a process* parameter values representing the effect of 
some physical, biological or chemical intervention, influence or event on the 
process. 

In the case of humans, the conditions to which the techniques are most ideally suited include 
conditions such as: 

Parkinson's Disease 

Schizophrenia 

Bipolar disorders / manic depression 
Cardiac disorders 
Myasthenia gravis 
Neuro muscular disorders 

Treatment of cancerous and tumorous cells and related disorders 
HIV / AIDS and other immune system disorders 
Hepatic disorders 
Athletic conditioning 
Treatment of pathogens 

Other disorders or diseases whose significant processes axe capable of being reduced to a 
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mathematical model 

An example of a mathematical model for Parkinson's disease is set out in Appendix B. In this 
regard, this is an example of a basic simple m lei, and is not intended to be limiting, but instead an 
5 example of the form of model that may be used. 

However, it will be appreciated that the process can be implemented with respect to any condition 
for which it is possible to construct a mathematical model of the condition. This is not therefore 
restricted to medical conditions, although the techniques are ideally suited for the application to 
10 conditions such as diseases or other medical disorders. 

The initial construction of a mathematical model, and subsequent determination of subject data in 
the form of parameter values may be achieved using any suitable technique, examples of which 
include: 

15 • Model Reference Adaptive Control (MRAC); 

• Neural Networks; 

• Complex Systems analysis; and, 

• Kalman filters. 

20 In this regard, the term Neural Networks will be understood to mean any method of creating 
mathematical models based on a supposed analogy with the functioning of the brain. In general, 
this is achieved by modelling the condition using a combination of non-linear functions including a 
suitably large number of terms to reduce mismatch errors. In this case, the model at a first 
approximation is a generic model that can be used to model almost any condition. This can then be 

25 trained through feedback, and associated restriction of parameters and the like, so that the model is 
constrained to model the respective condition of interest 

However, it will be appreciated that this is not intended to be limiting and any suitable technique 
could be used. 

30 

In any event, the system described above allows a mathematical model for the condition within a 
respective subject to be determined. Once this has been completed, the model is used to derive 
trajectories representing the progression of the condition under diflFerent circumstances, and in 
particular, under conditions of different medication or other treatment. 

35 

It will be appreciated that by calculating trajectories for the subject's current parameter and state 
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variable values, this allows the progression of the condition within the subject to be determined. In 
this instance, if the subject is already receiving medication, this allows the potential success of the 
existing medication regime to be assessed, or in the case of an untreated subject, the progression of 
the condition if the subject remains untreated. 

5 

However, in addition to this, it is also possible to determine and assess a number of different 
/ trajectories for a range of different parameter and state variable values. This allows undesirable 
/ trajectories to be identified and eliminated, for example because they do not lead to a stable end 
I point, or because they demonstrate chaotic behaviour, which would represmt rapid fluctuations in 
10, the condition of the subject. 

i 

i 

1 Once this has been con5>leted, it is possible to define stability sets, which correspond to sets of 
X.^ parameter and state variable values for which the remaining solution trajectories are acceptable, 
j Thus, this will correspond to trajectories that are non chaotic and result in the condition of the 
15 / subject reaching a suitable end point, such as a stable or periodic region, within which the condition 
I remains relatively xmchanged. 



20 



Having determined the stability sets, it is then possible to determine the changes required to the 
subject's current parameter and state variable values in order to ensure that the progression of the 
condition within the subject ends up following one of the acceptable trajectories. Thus, this is 
achieved by determining a control program representing the effects of external influences on the 
subject's parameter and state variable values. The control program is developed by mathematically 
modelling the external factors, such as the effect of medication, and this can be achieved using one 
or more of: 
25 • Liapunov Fimctions; 

• Dynamic Optimisation algorithms (such as ttie Euler-Lagrange Method); 

• Convex Set Algorithms (such as Kuhn-Tucker); 

• Any other suitable algorithm. 

30 The control program therefore represents a treatment program, such as a medication regime, which 
will cause the subject's parameter and state variable values to be modified towards those in the 
stability sets. 
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This therefore allows medication or the like to be prescribed on a case by case basis, to thereby 
ensure that subjects receive medication according to tailored regimes, thereby helping the subject 
progress to a condition state defined by one of the stability sets. 
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It will be appreciated that this will ideally lead to a stable end point or region representing the 
elimination of the condition. However, this is not always possible, and instead it may be necessary 
to select trajectories that '.^-id to a stable end point or region within which the symptoms and/or 
5 effects of the condition are eitl^er minimised, contained, or stabilised. 

Persons skilled in the art will appreciate that numerous variations and modifications will become 
apparent. All such variations and modifications which become apparent to persons skilled in the 
art, should be considered to fall within the spirit and scope that the invention broadly appearing 
10 before described. 
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Apendix A 

1 Problem statement: 

Consider an isolated coiitiguous colony of cells in the human body, of popu- 
lation x(t) cells for t > 0. This colony goes through successive generations 
of cells, whereby given sufficient resources, the rate of increase is proportion- 
ate to the existing population. However, there is competition for the finite 
resources (oxygen, blood etc.) provided by the colony's environment, and 
so rate of increase is increasingly retarded as the population expands. For 
a sufficiently large population, the rate of growth becomes negative, as cells 
die &om inadequate nourishment. 

This colony of cells secretes a substance S that is important for the well- 
being of the patient* Secreted by the colony into its external environment at 
a molar concentration this substance is eventually taken up by nearby 
protein transport mechanisms and used elsewhere. Rate of secretion takes 
the form of Michaehs-Menten kinetics, with the rate saturating to a maxi- 
mum once the colony's population becomes sufficiently large. Thus we can 
write a simple formidation of this system's dynamics as being 

x = x{t)iP-Rx{t))J{x{t)) (1) 

where P is the uninhibited population growth coefficient, R is the population 
retardation coefficient due to finite resources, Ki and K2 represent Michaelis- 
Menten quantities (more usually denoted Vm and K^f), /i is rate of uptake 
of S by transport mechanisms, and S {z) denotes the one-sided Heaviside 
function. Equation (1) is a more complicated form of the famous "logistic 
equation". 

2 Medical problem: 

A pathology is observed eunong a small percentage of the population, whereby 
y(f ) fluctuates sharply without apparent physical cause. In the medical com- 
munity this condition is known as "Y's Disease''. In case studies, it is ob- 
served that in some sufferers such fluctuations become progressively worse 
over a timespan of months or years, totally incapacitating the patient, until 
they become completely erratic, leading to death. 
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The three forms of medication at our disposal are 

• r. hormone (denoted ux) to stimulate the rate of unconstrfuned growth 
P of the cellulstr population; 

• a hormone (denoted uz) that strongly suppresses the rate of uncon- 
strained growth P of the cellular population (not conventionally used 
in treating Y's Disease); and 

• conventional medications to inhibit (denoted uz) or stimulate (denoted 
U4) the rate of uptake fi of y{t). 

Treatments involving these medications have not been particularly effective. 
Clinical trials using the growth hormone ui have, in some patients, unexpect- 
edly led to death, while control of ft has often been of fluctuating usefulness. 

3 Parameter identification, in vivo and in vitroz 

These dynamics are happeniiig within a living patient. With present tedmol- 
^Sy^ Pf ^1 and x(t) cannot be measured in vivo directly in any meaningful 
way; !€% can be estimated in tntro. y(t} can be measured non-invasively using 
fMRI, to a tolerance of dzSy. 

Instead of using conventional statistical methods of parameter estimation, 
which are virtually useless under the circumstances, we construct model equa- 
tions and a dynamic identifier algorithm (process illustrated in Figure 3A). 
Given a sufficiently large time-series of measurements {y (ti), y (*2) , - • - 1 y (*n)}i 
we are able to estimate the parameter values for oiu: specific patient to withhi 
an accuracy partly governed by the length of this dataset, as the model coef- 
ficients converge to those of the system. If no such convergence takes place, 
then this may indicate a structural disparity between modelling assumptions 
and the physical system. Diflferent forms of model equations can be em- 
ployed and cross-correlated (e-g, x ~ x(<) (P — i2 x^(t))S {x(t)), to confirm 
we have the best match of model to patient. Sufferers of the same medi- 
cal condition, at the same stage of the condition, should have models with 
a common structure, although the actual parameter values may vary firpm 
patient to patient. 

Having obtained these parameter estimates, they are substituted into 
the model equations to simulate the S3^tem. This is to chedc for validity, 
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by generating output from the model and correlating this with the ori^nal 
patient output. 

This identification process is rei>eated on a sufficiently large population 
of h -man or animal analogues to quantify the effects of the available medi- 
cations (listed above) on the system. As a consequence, it is found that the 
medicated system can be written 

X = x{t){{P + m^U2) -Rxit))J'(x(t)) (3) 

y = icf + ift)^ " (/^ + ^3 - U4) y(t)S{yit)) (4) 

where 

ttx € [0,0,15] (5) 

U2 € {[-l,-0.5],0} (6) 

tis 6 [0,0.2] (7) 

tt4 € [0,0.1]. (8) 

These constraints are imposed by drug safety concerns, or, in the case of i*2» 
the lack of fine control in delivering the hormone to the location. 

Ekiuations (3), (4) constitute the state equations of the system, more 
formally written 

z = f(z,u,A,t) (9) 
^here z = [r(*), y{t)r, u = [tii(t), «2(0, M*)^ Mt)^ and A = [P, Kx, K2, fi^. 

For our specific patient, we find the reference model of the identifier 
algorithm converges asymptotically to ODEs of the form of equations (1), 
(2), with normalized parameter values R = 1.6; P = 2; Kx = 10; K2 = 0.8 
and /X = 1.2 when the patient is immedicated, i.e. u = O. 

4 Global stability analysis 

As illustrated in Figure 3B, a global stability analysis is conducted on the 
structure of equations (1) amd (2), across a range of values for P and R 
(Figure 11), The top illustration of Figure 11 is a convergence map, of 
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repeated iterations on x(0)=l. Dark grey denotes stability (the sequence 
converges eventually); black denotes extreme stability (the sequence tends 
immediately to a final value, usually extinction) and white denotes no final 
convergence after a designated number of iterations - in this case, five thou- 
sand - either due to the sequence diverging, or else the onset of chaos. The 
bottom illustration is a false-colour contrast enhancement of this convergence 
map, rendered here in B&W. 

To measure the dependence (if any) of these maps on our initial choice 
of x(0), a second, "chequerboard'' set of convergence maps is generated, of 
x(0) against P for fixed values of it This is depicted in Figure 12, 

It can be seen from the structure of Figm-e 11 that 

1. This system is susceptible both to chaos and to extinction; 

2. That, consequently, the manifestation of strong fluctuations in y(t) in 
**Y*s Disease" need not necessarily be a disease at all, but a peturba- 
tion of the dynamics of an otherwise healthy system into a region of 
instability or extinction; 

3- That, therefore the unexpected deaths of patients under ui growth- 
hormone therapy need not imply any inherent toxicity of the hormone, 
but be due to the inducing of chaos or even the extinction of the system 
evident for P + ui being pushed to too high a value (e.g, P + tii > 3); 

4. That, given the complex population dynamics of the underlying ceU 
colony, attempts to control y(t) without controlling x(t) are futile; 
hence, modification of the cell population using hormones txi, U2 is fun- 
damental to medicating the patient. Contrary to expectation, conven- 
tional medications and U4 that alter uptake of y(t) can only ever be 
of limited practical use in regulating y{t); 

5. The most stable region of this system is in neighbourhoods surrounding 
{P = 1,R <^ 2}, although e'ven liere there exist dangerous zones of 
chaotic structure (marked white in upper illustration and black in lower. 
Figure 11). 
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For specified values of P and 12, there exists two pairs of equilibrium 
points [a;*,2/*]'^ such that 

i{x\y^) = y{x\y^)^Q, (10) 

namely, 

KyiT = (o,o). (11) 

and 

The first of these points corresponds with the trivial equilibrium of extinction; 
the second, of a dynamic non-zero equilibrium of the population (visible in 
the frames of Figure 12 as a line of black dots). 



5 Single-Liapunov control 

The next task is to control the system, steering it into a stable configuration 
to improve the patient^s condition (the process depicted in Figures 4A 
and 4B). An early step of this is to decide upon a desired target set ri. 
As illustrated at various levels of magnification in Figure 13, this partly 
depends upon the scale of the system. As well as the large-scale pool of 
stability visible for P € [0.8, 1.19], even within the diaotic domain P € 
[1.2, 1.4] small locally-stable neighbourhoods appear to exist. 

In the case of this patient, the cell colony is to have its growth param- 
eter permanently translated to P 1, by setting uz{t)\ = — 1 Vf. In the 
XY plane for P = 1, = 1.6, we then define ti to be a neighbourhood 
of the stable equilibrium point [:>:|i2/2]^* This is illustrated in Figure 14, 
where the relative stabiUiy of points [a;(*), y(*)]''^ is mapped left in the usual 
grey-scale (black = extreme convergence, white = failure to converge) and 
right using false-colour shading to enhance difierences in relative stabihty, 
rendered here in B&W. The equilibrium point [xj, y^]"^ is tagged by a dark 
point, and can be seen to lie on a thin curve of relatively stable points. The 
elliptical neighbourhood that constitutes ri is shown. All points that reach 
this neighbourhood will converge to [a:^, y^]^- 

A Liapunov function 14 is then constructed around the geometry of Ti, 
and a control program computed using VVi-f(z,u, A, t) < 0, to use the 
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constrained ux, U3 and U4 to steer the patient's state to Tx. Candidate so- 
lution trajectories are computed and presented to the clinician (e.g. see 
Figure 15a, for a solution trajectory issuing from the point [x(0),i/(0)]'^ = 
[0.001,0]'^). If these solution trajectories or p* are in some sense unsatis- 
factory, quantitative optimization is performed (e.g. Figure 15b shows a 
slightly different trajectory steered fr-om the same initial conditions to the 
target set, in which some robustness is sacrificed in order to abandon any 
reliance upon the medications t/3, 1*4). 

The set of all z(0) from which such Vi-derived control is possible is comr 
puted and denoted R?. 

6 Double-Liapunov control 

Now consider the problem where there exists tmcertainty about the precise 
value for R. This may be due to background fluctuations preventing the 
identification algorithm from attaining a more precise estimate of R than 
R±SR, or it may be that the patient's biology engages in a potentially adverse 
response when P is modified via hormone therapy, whereby R actually is 
shifted by amount SR, Perhaps this is due to the disease, or possibly it is 
a natural body mechanism. Either way, it is necessary to repeat the control 
process of the previous section in such a way that it is guaranteed to be 
robust against such fluctuations (the process depicted in Figures 6, 7a and 
7b). 

The equations (3), (4) are re-written, 

X = x(t)((P + m-U2) -(iJ + «;)x(*))J~(x(t)) (13) 

if = T^?%-rH*))-(/* + «3-u4)y(0J"(y(0) (14) 

where w G W^] and w'^ denotes the maximum possible fluctuation in 

R, estimated for our human patient to be = 0.1. 

In this "Game against Natiu'e'' it is assumed that these fluctuations are 
deployed intelligently and aggressively against our control strategies by a 
player designated **Nature'', in an attempt to disrupt the system- The task is 
then to design control strategies that are robust against any such interference. 
The target set and Liapunov function Vi^ are taken to be the same as 
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Tu Vi' To decide upon appropriate choices for ra. Figure 11 and Figure 
12 are reexamined. Even at R ^ 1.1, it is clear that for P = 1. domains 
corresponding to crtinction cannot easUy be readied by Nature directly; 
the most inunediate way to disrupt population dynamics is by driving x(t) 
either to x(t) 0 or x(t) -> :> 1.2. A value of = 1.4 is chosen, and a 
neighbourhood around this nominated as r^. 

An appropriate Liapunov function V2 is designed, and Nature strategies 
P- generated such that VVi.f < 0 and p** attempts to obstruct our con- 
trol efforts. The control strategy p— is computed, based on the criterion 
yVx^-f < 0. The sets Wt (points [x(0),y(0)f from whidi p— steers tra- 
jectones safely into r,„. against aU possible ui) and W2 (points [x(0), y(0)J'^ 
from which p** steers trajectories into despite all possible u) are plotted. 
As before, given lx(0), y(0)]'', candidate solution trajectories are computed 
under p— and presented to the dLnician (e.g. see Figure 16. for a sohition 
teajectory issuing fi^m the point (x(0), y(0)P = [1. 0]-*-. giving successful con- 
trol programs for u despite the counterattack 11/). These control 
can be subjected to another layer of optimization, or else the constraints 
[u , u-*-], [to , w+J modified to- improve clinical control of the patient. 

Thus it has been demonstrated that these tools can be emplqyed to trans- 
form the medication <rf patients through (a)dynamic parameter and state 
estimation. (b)the calculation of single-function control and (c) the compu- 
tation of robust control in the presence of random or activriy hostile elements. 
In most cases such medication control strategies will be profoundly^ different 
from those of conventional ^npirical medication. 
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Appendix B 

The equations set out below represent a simple 3-equation model of exogenous dopamine 
pharmacokinetics in the synaptic cleft of a dopaminergic neuron. 

5 



10 



where: 

[^cLa levodopa concentration in bloodstream 

ku ' ' l^v^^^P^ transport parameters across blood-brain barrier 

15 //^^ : "sink" coefficient for levodopa in bloodstream 

[^elyn • l^vodopa conccntration in synaptic cleft 
kiT' » k*iT' ' conversion parameters for dopa decarboxylase 
Mesyn • "sink" coefficient for levodopa in cleft 
\f^e lyn • dopamine concentrarion in synaptic cleft 

20 , : dopamine transport coefficients for neuron mass 

Mdsyn ' "sink" coefficient for dopamine in cleft 
{d^): mean concentration of signalling dopamine 
AyWy^: control variables for introducing levodopa into blood 
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Accordingly, it will be appreciated by persons skilled in the art that the above mentioned equations 
can be used in determining a model representing the progression of Parkinson's disease within a 
subject. 



